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Abstract. We present three 14400 square degree relative extinction maps of the Galactic Plane (\b\ <20°) obtained from 
2MASS using accumulative star counts (Wolf diagrams). This method is independent of the colour of the stars and the variation 
of extinction with wavelength. Stars were counted in 3.^5 x 3.^5 boxes, every 20'\ 1° x 1° surrounding fields were chosen for 
reference, hence the maps represent local extinction enhancements and ignore any contribution from the ISM or very large 
clouds. Data reduction was performed on a Beowulf-type cluster (in approximately 120 hours). Such a cluster is ideal for this 
type of work as areas of the sky can be independently processed in parallel. We studied how extinction depends on wavelength 
in all of the high extinction regions detected and within selected dark clouds. On average a power law opacity index (f3) of 1.0 
to 1.8 in the NIR was deduced. The index however differed significantly from region to region and even within individual dark 
clouds. That said, generally it was found to be constant, or to increase, with wavelength within a particular region. 
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' 1. Introduction 

Dust is not only one of the most important coolants in molec- 
ular clouds but is also an excellent tracer of molecular hydro- 
^ gen. Determining its distribution is thus crucial to understand 
■ the early stages of star formation. It gives us information, for 
" example, on how clouds fragment and ultimately this knowl- 
edge must be reconciled with the initial mass function (Padoan, 
Nordlund, & Jones ll997t . Moreover by combining dust extinc- 
tion maps with infrared and submillimeter dust emission maps, 
we can derive basic knowledge of grain properties and search 
for evidence of grain evolution. 

Extinction towards an individual star can of course by de- 
termined by measuring its colour excess although this method 
requires knowledge of the star's intrinsic colour (e.g. He et al. 
119951 Racca et al. 12002ft . The technique is also highly selec- 
tive as it only gives us information along very precise lines 
of sight. An alternative approach, developed many years ago 
(Wolf 1 19231 Bok ll956V and particularly appropriate for dark 
clouds, is to count the number of stars in the cloud's vicinity 
down to a limiting magnitude and compare this with a control 
"unextincted" region nearby. Such so-called Wolf Diagrams are 
of limited use in the optical, as one can only probe the outer 
peripheries of a cloud where the extinction is low. The method 
works well however in the infrared, due to its improved trans- 
mission with respect to the optical, allowing regions with Ay 
values as high as 30-40 to be measured. 
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Dust not only absorbs light from background stars but also 
reddens it. By measuring extinction in a number of wavebands, 
one can determine how it varies with wavelength (A). Normally 
this variation is assumed to be a power law, i.e. Aa cx 
where P is the opacity index. A constant value for p however 
is not necessarily guaranteed a priori. 

One drawback in the past of the star counting method, at 
least for large areas, is that it is computationally very expen- 
sive. This type of data reduction, however, is ideally suited 
to a parallel computing environment provided, for example, 
by a Beowulf-type cluster. Large swathes of the sky can be 
processed simultaneously and independently of each other. 
Moreover, the possibility for exploiting such "number crunch- 
ing" facilities coincides with the availability of large ground 
based IR surveys such as 2MASS and DENIS. For the first time 
its possible to manufacture extinction maps covering large frac- 
tions of the sky. 

Here we present a high spatial resolution relative extinction 
map of the Galactic Plane, in J, H, and K, derived from the 
2MASS database using the star count method. The method and 
its limitations are described in detail in Sect. 2, the determina- 
tion of the noise is discussed in Sect. 3. Results of the opacity 
index computations are presented in Sect. 4 and our conclu- 
sions are given in Sect. 5. 

2. Method and Limitations 

Relative extinction maps were determined using Wolf diagrams 
(e.g. Kiss et al. l2000ll . At each position accumulated star counts 
were performed in a small box. These counts were compared 
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Fig. 1. JHK colour composite of the star count map (top), relative extinction map obtained from the J-band data (middle) and 
three sigma noise due to the non-uniform distribution of stars in the relative extinction maps for J (bottom). The noise is displayed 
in linear scale from zero (white) to 0.7 mag (black) of optical extinction. The rectangle marks the region magnified in Fig.|2l 



to the average accumulated star count in a (co-centred but 
larger) comparison field to determine the relative extinction. 
The method is independent of stellar colour and variation of 
extinction with wavelength. A number of assumptions, how- 
ever, have to be made (e.g. Froebrich et al. 12005b : (1) stars are 
distributed uniformly and all apparent voids are due to extinc- 
tion (2) one may define a constant average absolute magnitude 
for the stars in a box (3) the completeness limit of the catalogue 
does not depend on position. In reality these assumptions may 
only be valid under certain conditions. For example reference 
to an average absolute magnitude is only meaningful if the box 
contains enough stars. Moreover stars are not distributed uni- 
formly on the sky but if the control field is not too far away, the 
assumption of uniformity is locally valid. Thus, to apply the 
method properly one must ensure that the choice of box size, 
control field, etc., are within certain limits (see below). 

Using the 2MASS catalogue already ensures a high level of 
uniformity and quality of the photometry. To avoid photometric 
errors for faint objects we only selected sources with a signal to 
noise (S/N) ratio > 5. This corresponds to quality flags A/B/C 
in the 2MASS catalogue, and an error in stellar magnitude of at 
most 0. 17 mag. Such a value is slightly higher than our stepsize 
in brightness (counts were performed every 0. 1 mag). However, 
the extinction is derived using an average for all stars in the box 
and is thus determined with greater precision. 

As we have emphasised, accumulated star counts are only 
meaningful, providing boxes contain a sufficient number of 
stars. In particular the S/N in the final relative extinction maps 
strongly depends on the number of stars in the area being 
counted and hence on the box-size of the reseau (see Appendix 
C in Froebrich et al. 12005 I f or a detailed discussion). Taking into 
account the density of stars in the 2MASS catalogue in all 3 fil- 
ters, we selected a box-size of 3.'5 x 3.'5 to perform the counts. 
This ensures approximately 25 stars per box (S/N = 5, assuming 
a Poisson distribution) at the completeness limit for 99.5, 73.5, 



63.5 % of the area with |6| < 10°, in J, H, and K, respectively. It 
also leads to a mean density of about 100 stars or more per box 
in the J-band for a large area (^ 2800 0°). Such a density con- 
verts to around one star per 20''x20''. Hence performing star 
counts with a spatial frequency of 20'' ensures optimal sam- 
pling of the 2MASS catalogue. Note that this oversampling of 
about 10 is not unconditionally required in regions with fewer 
stars, but was kept to ensure a uniform pixel size throughout 
the whole extinction map. 

Counts in the 3. '5 x 3. '5 box where compared with a co- 
centred larger control field. For the latter we chose a l°xl° 
area. This has the effect that our maps record local extinction 
enhancements, due to small clouds or globules, and not global 
extinction effects due to the ISM or extended dark clouds. 

Achieving a high S/N in our maps depends on extending 
the star counts to the faintest possible magnitude (and thus the 
largest possible number of stars). Care however has to be taken 
as the completeness limit of the catalogue varies strongly with 
galactic coordinates. We obtained this limit for each position by 
determining the peak in the histogram of stars per magnitude 
bin in the associated control field. Our relative extinction maps 
are then obtained exclusively using stars 0.3 mag brighter than 
this completeness limit. 

The exercise of counting stars can be described as an em- 
barrassingly parallel problem, as results from one part of the 
sky do not depend on any other. The entire 2MASS catalogue 
covering \b\ < 20° from the Galactic Plane was processed with 
a sampling frequency (effective pixel size) of 20'' using our 32- 
node Beowulf-type cluster of 2.67 GHz P4 HyperThreat pro- 
cessors. The three 14400n° relative extinction maps (in J, H 
and K) took about 40 hours each to calculate. 

In Fig. [l] we present a three colour composite (J, H and K) 
full size star count map (top panel) and the relative extinction 
map (middle panel) obtained from the 2MASS J-band data. 
The small rectangle in the figure marks the region blown-up 



D. Froebrich et al.: A Galactic Plane Relative Extinction Map from 2MASS 3 




303 302 301 300 303 302 301 300 303 302 301 300 

m 



Fig. 2. Relative extinction map obtained from J-band 2MASS data of the Musca Dark Cloud. Contours in the left panel start 
at Av = 1 mag and are in steps of 1 mag. We marked previously known and new regions with peak extinction values larger than 
Av = 2 mag. "Fake" globules due to bright stars are labelled with *, and numbers represent regions from Vilas-Boas et al. d 199411 
(MU 16 could not be detected). The middle panel shows the relative extinction map in grey-scale and the right panel the 100 /im 
IRAS image of this region. In the marked rectangle, around NGC 4372, extinction values are less reliable due to additional noise 
(see also Sect.|2|). 



in Fig.|2| The latter figure, of the Musca Dark Cloud, demon- 
strates the full resolution obtained in our relative extinction 
maps (in contours in the left and grey-scale in the middle 
panel). Known dark clouds are labelled, as well as "fake" glob- 
ules (see below). The panel to the right of Fig.|2|shows for com- 
parison the IRAS 100 /im image and nicely demonstrates the 
correspondence between dust in emission and absorption. 

Comparing accumulated star counts in a 3.^5 x 3f5 box with 
the surrounding 1 ° x 1 ° field is ideal for finding small regions 
of enhanced extinction. All fields containing high extinction 
regions, which are more extended than a significant portion of 
one square degree, however, will not have the correct extinction 
values. In particular this includes the large Orion, Taurus, and 
Ophiuchus clouds. Here a larger comparison field is needed, 
but see the discussion in Sect.|4| 

Very rich star clusters (mainly globular clusters) are prob- 
lematic in the sense that star counts are anomalously high at 
their centre compared to their peripheries. This gives rise to 
regions of apparent negative extinction at the cluster centre 
and apparent positive extinction further out. Fig.|2l showing the 
globular cluster NGC 4372, illustrates this effect. Fortunately 
these regions are usually not very close to areas of real higher 
extinction. Another effect is that extremely bright stars prevent 
detection of sources close by and hence also mimic high extinc- 
tion regions (also see Fig.|2|for a few examples). Finally there 
are physical and effective gaps in the 2MASS catalogue which 
also result in higher apparent extinction. Fortunately these re- 
gions cover only 0.006 % of the whole sky. 



3. Noise Determination 

We determined the noise in our images to compare our method 
with the NICER technique of Lombardi & Alves SZOOTi . Our 
relative extinction maps show noise due to the non-uniform dis- 
tribution of stars, a I , and the box- size within which we count 
stars (so-called dither noise), (72. The latter manifests itself as 
structures of around the size of the box used for star-counts. 
We thus smoothed the original maps with a filter of width the 
box- size and subtracted the resulting image from the original 
to obtain the noise due to the non-uniform distribution of stars. 
These noise images were then gaussian filtered (5' FWHM), 
in order to obtain the same resolution as in Lombardi & Alves 
(1200111 . and the 3a noise level was determined. The noise is 
transformed into Ay using the conversion factors 4.04, 6.47, 
9.75 for J, H, and K, respectively, given in Mathis (1 199011 . The 
bottom panel in Fig.Qshows in grayscale the resultant ai noise 
in our J-band map, averaged over 1° x 1°. 

This noise ranges from 0.2- 1.0 mag (J), 0.4-1.8 mag (H), 
and 0.6-3.8 mag (K) in equivalent optical extinction. The val- 
ues for the J-band data are slightly higher than obtained by the 
NICER technique (compared with the maps of Lombardi & 
Alves (l2QQni in the Orion region at about / = 210°, 6= -18°). 
All three noise maps show the same principle structure 
with the lowest noise levels obtained north and south of 
the Galactic Plane within ± 60° longitude from the Galactic 
Centre. Regions directly in the Galactic Plane and within the 
large dark clouds (Ophiuchus, Orion, Taurus) are affected by 
higher noise, due to extended extinction and hence a smaller 
number of stars. 

The additional noise, due to dithering, is only slightly de- 
pendent on the chosen box-size. A larger number of stars in the 
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box reduces <J2 . By varying the box-size the value for this noise 
can be estimated. For our box-size we obtain typical values of 
about 0.7, 1.0, 1.3 mag optical extinction for the J, H, K-band 
dither-noise, respectively. See the Appendix of Froebrich et al. 
(12005ft for a detailed description of this effect in the context of 
a particular region (IC 1396). The actual detection limit in our 
relative extinction map can be determined by adet = v^cr^ + cr^ 
for each position and filter individually. 

4. NIR Opacity Index 

Using our three large scale maps we can determine how extinc- 
tion varies in the NIR. The power law index (3 of the wave- 
length dependence of the extinction can be determined using: 

(3x,x, = log(AA,/AA,)/log(A2/Ai) (1) 

where Aa^ is the extinction at a wavelength Ai. We chose 
the wavelength for the 2MASS filters^ of 1.235, 1.662, and 
2.159 /im for J, H, and K, respectively. Since we have three fil- 
ters there are three possible combinations to determine (3 (JH, 
JK, HK), two of which are independent. 

Care has to be taken in deriving /3, since it is very sensitive 
to small errors in extinction especially at low extinction values. 
Hence only regions which are three sigma above the local noise 
level were considered. A further more complex source of error 
is the size and position of our control field. Due to extinction, 
a fraction / of the control field might just show a fraction e of 
stars, compared to unextincted regions. This leads to a mean 
extinction Ac in the control field and hence a smaller apparent 
measured extinction AA=Areai— Ac- It can be shown easily 
that under the assumptions / • (1— e) <C 1 and Ac <^ Areai, the 
derived /3 value is not influenced, even if the measured extinc- 
tion values are wrong by Aq. The size of Ac is smaller than 
the noise in our maps as long as a/TV • / • (1— e) < 1, were TV 
is the number of stars in our small box. According to our box- 
size and the star density, this holds for most regions as long as 
/ is smaller than 0.3. Considering this we also excluded all re- 
gions where more than 30 % of the control field is influenced 
by extinction larger than the three sigma noise from the opacity 
index determination. 

The distribution of /?jh, determined over the whole area 
of our map is rather broad and peaks between 1.0 and 1.8. 
The same wide distribution is found for /^hk, but at slightly 
larger P values (shifted by +0.1). The reason for this particu- 
larly broad distribution can be found when we look at groups of 
dark clouds in detail. We find they can be categorised accord- 
ing to whether the opacity index 1) is roughly constant, 2) in- 
creases, or 3) decreases with wavelength in the NIR. Note that 
even within individual small clouds (3 can vary significantly. In 
Table[l]we list peak positions and width of the distribution of (3 
for a selection of dark clouds. Most regions belong to Category 
1 or 2, although a small number (e.g Cepheus) can be classified 
as Category 3. Note that (3 varies significantly and ranges from 
1.0 (e.g. S 140) to 2.0 (e.g. Circinus). 

^ http://www.ipac.caltech.edu/2mass/releases/allsky/doc/explsup.hti 



Table 1. Peak position and width of the distribution of the 
opacity index (3 for selected dark clouds, determined between 
J and H (/3jh) and between H and K (^^hk)- In some regions 
/?HK could not be determined properly, due to low S/N. 



Region 






Region 




PuK 


CrA 


L4±a3 


1.4±0.4 


LupusIII 


L5zb0.3 


L9±0.4 


Ophiuchus 


1.5zba3 


L7±0.4 


LupusIV 


L5zba3 


1.9dba3 


LupusI 


L5±a3 


L7±0.5 


Monocerus 


L4±0.5 


L0±0.4 


LupusII 


1.3di0.2 




Cepheus 


L4±0.4 


L2±0.4 


Musca 


1.7±a3 




S140 


L0±a3 


L5±0.4 


Circinus 


L7±0.4 


2.0±0.4 


Serpens 


L7±0.4 


2.0±a3 



5. Conclusions 

We have shown that accumulated star counts in the NIR can 
be used to obtain large scale relative extinction maps. Parallel 
techniques allowed us to process 14400 square degrees of the 
sky within 40 hours computing time on a 32-node Beowulf- 
type cluster of 2.67 GHz P4 HyperThreat processors. The accu- 
mulated star counts are performed independent of stellar colour 
and only local extinction enhancements are determined. Hence 
the wavelength dependence of the dust extinction within the 
dark clouds can be investigated. A study of selected clouds 
shows that the opacity index P ranges from 1-2. In most cases a 
constant or increasing opacity index with wavelength is found. 
The large scatter of P, even within a particular cloud, leads 
to the conclusion that a uniform opacity index should be used 
with great care for extinction corrections within dark clouds. 
This is in contrast to the general ISM where mixing of the dust, 
plus averaging over long lines of sight, ensures a single index 
suffices. 
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